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Abstract 

We present measurements of pion elliptic flow (V2) in Au+Au collisions at -\Js NN = 200, 62.4, 39, 
27 and 19.6 GeV, as a function of event-by-event charge asymmetry (A ± ), based on data from the 
STAR experiment at RHIC. We find that ii~ elliptic flow linearly increases (decreases) with 
charge asymmetry for most centrality bins and for all the beam energies under study. The slope 
parameter (r) from V2(A±) difference between n~ and n + shows a centrality dependency similar to 
calculations of the Chiral Magnetic Wave. The measurements of charge separation with respect 
to the reaction plane in search of Local Parity Violation and the Chiral Magnetic Effect are also 
presented for Au+Au collisions at yfs^ = 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, and for 
U+U collisions at 193 GeV. 



1. Introduction 

In heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC), energetic spectator protons produce a strong magnetic field peaking around 
eB y « m% U. The interplay between the magnetic field and the quark-gluon matter created in 
the collisions is characterized by two phenomena: the Chiral Magnetic Effect (CME) and the 
Chiral Separation Effect (CSE). The CME is the phenomenon of electric charge separation along 
the axis of the magnetic field in the presence of a finite axial chemical potential (e.g. chiral po- 
tential due to fluctuating topological charge) dHHSH. STAR HQ] and PHENIX HH col- 
laborations at RHIC have reported experimental observations of charge asymmetry fluctuations 
possibly providing an evidence for CME. This interpretation is still under intense discussion, see 
e.g. llU lllll and references therein. The Chiral Separation Effect (CSE) refers to the separation 



of chiral charge along the axis of the magnetic field at finite density of vector charge (e.g. electric 



charge) \YM3- 

In a chirally symmetric phase, the CME and CSE effects form a collective excitation, Chi- 
ral Magnetic Wave (CMW), a long wavelength hydrodynamic mode of chiral charge densi- 
ties 11141 llal . The CMW manifests itself in a finite electric quadrupole moment of the colli- 
sion system, where the "poles" ("equator") of the produced fireball acquire additional positive 
(negative) charge [14]. This effect, if realized, will be reflected in the measurement of charge- 
dependent elliptic flow. Elliptic flow, characterized by a second-order harmonic in the particle 



A list of members of the STAR Collaboration and acknowledgements can be found at the end of this issue. 
Preprint submitted to Nuclear Physics A November 16, 2012 



azimuthal distribution (0) and quantified by the Fourier coefficient V2 lfl6tl . refers to the collective 
motion of particles with respect to the reaction plane (i^rp): 



v 2 = <cos[2(0-^rp)]>. (1) 



Taking pions as an example, on top of the baseline vS* 86 ^*), a CMW will lead to M1411 



V2(**) = + (^)A ± , (2) 

Pe 

where q e , p e and A ± = (N+ - N_)/(N + + iV_) are the quadrupole moment, the net charge density 
and the charge asymmetry of the collision system, respectively. As (A ± ) is always positive, 
A ± -integrated V2 of n~ should be above (below) the baseline due to CMW. However, the 
baseline V2 may be different for n + and ii~ in the first place because of several other possible 
physics mechanisms 017lll8lll9ll . so it is less ambiguous to study CMW via the A ± dependency 
of pion v'2 than A ± -integrated v>2. 

In Sec. 2, we present A ± -differential measurements of pion \>2 for Au+Au collisions at 
y/SffN =200, 62.4, 39, 27 and 19.6 GeV. We find that pion v'2 exhibits a linear dependence on A ± , 
with positive (negative) slopes for if~ (tt + ). The slope difference between if~ and 7r + is studied as 
a function of collision centrality. In Sec. 3, the measurements of charge separation correlator are 
used to search for CME and Local Parity Violating (LPV) effects. 



2. Pionv 2 (A ± ) 




Figure 1: (Color online) The example of 30-40% Au+Au collisions at 200 GeV t22ll . (Left) Pion v'2 { 2 1 as a function 
of observed charge asymmetry. (Right) v>2 difference between n~ and n + as a function of charge asymmetry with the 
tracking efficiency correction. The errors are statistical only. 

Charged particle tracks were reconstructed in STAR TPC [20], with pseudorapidity cut \rf\ < 
1. The centrality definition and track quality cuts are the same as in Ref. 112 111 , unless otherwise 
specified. This study is based on Au+Au samples of 200M events at 200 GeV from RHIC year 
2010, 60M at 62.4 GeV (2010), 100M at 39 GeV (2010), 40M at 27 GeV (2011) and 20M at 
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19.6 GeV (201 1). All were obtained with a minimum-bias trigger. Only events within 40 cm of 
the center of the detector were selected for this analysis. In the calculation of charge asymmetry, 
(anti)protons with transverse momentum pj < 0.4 GeV/c were excluded to reject beam pipe 
protons. A distance of the closest approach (dca) cut (< 1 cm) was also applied to reduce the 
number of weak decay tracks or secondary interactions. To select pions, we eliminate charged 
particles 2<x away from the expected TPC energy loss for pions. 

Elliptic flow measurements were carried out with the two-particle cumulant method V2{2} lfl6l 
I2H1 for 200 and 62.4 GeV, and \'2{i] sub) approach for the rest beam energies, where two subevents 
consist of charged particles with 77 > 0.3 and 77 < -0.3, respectively. Pions at positive (negative) 
77 are then correlated with the subevent at negative (positive) 77 to calculate \>z- The 77 gap of 0.3 
unit suppresses short-range correlations such as Bose-Einstein interference and Coulomb final- 
state interactions lElll . The 77 gap was also used in the V2{2) analysis in a similar way. To focus 
on the soft physics regime, only pions with 0.15 < pr < 0.5 GeV/c were used to calculate the 
/?7-integrated V2. Taking 30-40% 200 GeV Au+Au as an example 12211 . we show pion V2 as 
a function of observed charge asymmetry in the left panel of Fig. Q] n~ V2 increases with the 
observed A ± while 7r + V2 decreases with a similar magnitude of the slope. After the tracking 
efficiency correction for the charge asymmetry, the vi difference between n~ and ti + is fit with 
a straight line in the right panel. The slope parameter r, or 2q e /p e from Eq. 2, is positive and 
qualitatively consistent with the expectation of the CMW picture. The intercept of the linear fit 
is non-zero, indicating the baseline v'2 for if and n + are different. 
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Figure 2: (Color online) The slope parameter r, supposedly 2q e /p e , as a function of centrality. For comparison, we also 
show the UrQMD [23] simulation for 200 GeV Au+Au, and the calculations with CMW 1 14] with different duration 
times. The grey band represents the systematic uncertainty due to varied dca cuts and the tracking efficiency. 

In both V2{2) and V2{?7 sub), there are correlations not related to the reaction plane, and not 
suppressed by the 77 gap, for example due to back-to-back jets. They are largely canceled out in 
the V2 difference between n~ and n + . Correlations between daughters of weak decays like A/A 
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may not be canceled, and may contribute to the intercept of Vz(A ± ) difference. However, for 
neutral particles like A, the decay daughters on average will not contribute to the numerator of 
observed A ± , so they will not create a correlation between V2 difference and A ± out of nothing. In 
other words, this effect will not change the slope parameter from zero to finite. The denominator 
of observed A ± will be increased due to this effect, and thus the observed slope parameter will be 
increased by a scale factor, related to the A/x ratio, which requires further systematic study. 

We follow the same procedure as above to retrieve the slope parameter r for all centrality bins 
and all the collision systems under study. The results are shown in Fig. [2] together with the sim- 
ulation calculations with UrQMD event generator [23] and with the theoretical calculations with 
the CMW effect with different duration times for the magnetic field lfl4ll . For most data points, 
the slopes are positive and reach a maximum in mid-celtral/mid-peripheral collisions. The slopes 
extracted from UrQMD events are consistent with zero for 15-60% collisions, where the signal 
from the real data is prominent. On the other hand, the CMW calculations demonstrate a similar 
centrality dependency of the slope parameter, though quantitative comparison between data and 
theory requires further works on both sides to match the kinematic regions of the analyses. 



3. Charge separation with respect to the reaction plane 
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Figure 3: (Color online) The three-point correlator, y, as a function of centrality for Au+Au collisions from 200 GeV to 
7.7 GeV [24]. For comparison, we also show the results for Pb+Pb collisions at 2.76 TeV l25ll . The errors are statistical 
only. 

The concept of Local Parity (V) Violation (LPV) in high-energy heavy ion collisions was 



brought up by Lee et al. [26, 27, 28] and elaborated by Kharzeev et al. B29I1 . In non-central 
collisions such a !P-odd domain can manifest itself via preferential same charge particle emission 
for particles moving along the system's angular momentum, due to the Chiral Magnetic EffectO, 
3]. To study this effect, a three -point mixed harmonics azimuthal correlator was proposed 13011 : 



y = <cos(0„ + (pp- 2iArp)>, 
4 



(3) 



where a and /3 denote the particle type: a, y3 = +, -. The observable y is !P-even, but sensitive 
to the fluctuation of charge separation. STAR measurements of the correlator were reported 
for Au+Au and Cu+Cu collisions at 200 GeV and 62.4 GeV Jf^O], showing the clear difference 
between the opposite sign and the same sign correlations, qualitatively consistent with the picture 
of CME and LPV. Fig.[3]presents the extension of the analysis to lower beam energies at RHIC. 
The STAR results are based on Au+Au samples of 57M events at 200 GeV from RHIC year 
2007 Q, 7M at 62.4 GeV (2005), 100M at 39 GeV (2010), 40M at 27 GeV (2011), 20M at 
19.6 GeV (2011), 10M at 11.5 GeV (2010) and 4M at 7.7 GeV (2010). For comparison, we also 
show the results for Pb+Pb collisions at 2.76 TeV 12511 . A striking similarity exists between 200 
GeV Au+Au and 2.76 TeV Pb+Pb, and a smooth transition occurs from 200 GeV to lower beam 
energies starting from the peripheral collisions. 
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Figure 4: (Color online) The difference between the opposite sign and the same sign correlations as a function of central- 
ity for Au+Au collisions from 200 GeV to 7.7 GeV 1 24] . For comparison, we also show the results for Pb+Pb collisions 
at 2.76 TeV [25]. The errors are statistical only. 

Initially it was expected that the opposite sign (yos) and the same sign (yss) correlations 
would be symmetric around zero due to the charge separation induced by LPV and CME. How- 
ever, there could be common physics backgrounds in both correlations. For example, in central 
collisions the strong radial flow tends to push particles to the same direction regardless of the 
charge sign, and that effect will reduce both 70s and yss by the same amount. In peripheral 
collisions, the multiplicity is smaller and the system is more influenced by momentum conserva- 
tion, which tends to increase both correlations in the same way. Also, the statistical fluctuation 
of the correlator could be larger out-of-plane than in-plane due to the geometry of the collision 
system, which contributes a negative background. To reduce such mutual backgrounds, we take 
the difference between yos and yss as the signal, shown in Fig. [4] The signal persists almost 
unchanged up to 2.76 TeV and down to around 11.5 GeV, and seems to disappear at 7.7 GeV. 
To be more conclusive on the transition of the signal, more statistics are needed for collisions at 
11.5 and 7.7 GeV. 

One major background in yos -yss comes from processes in which particles a and (3 are prod- 
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Figure 5: (Color online) Comparison of multiplicity (left) and pj -integrated vi (right) between 200 GeV Au+Au and 
193 GeV U+U. The errors are statistical only. 



ucts of a cluster (e.g. resonance, jet, di-jets) decay, and the cluster itself exhibits elliptic flow (V2) 
or different emission between in-plane and out-of-plane [31]. The corresponding contribution to 
the correlator can be estimated as 0211 

^clust/event-^pairs/clust . _ . .. 

(cos((p a + 4>p- 2^ dust )) c i ust • v 2 ,ciust, (4) 

*V pairs/event 

where (...) c i us t indicates an average over pairs consisting of two daughters from the same cluster. 
In the picture of global momentum conservation [33] and/or charge conservation D~4Tl . the whole 
system may be considered as one cluster. To estimate this type of V2-related background, we 
proposed to study collisions with sizeable V2 and almost no magnetic field via very central U+U 
collisions Bill . The U+U collisions at 193 GeV were carried out in RHIC year 2012, and the 
data taken by STAR contain both the minimum-bias trigger and a dedicated online trigger to 
take events with 0-1% spectator neutrons, so that the very central U+U collisions were selected 
with the magnetic field greatly suppressed. Uranium is larger than gold, which is reflected in 
the multiplicity distribution in the left panel of Fig. \'2{t] sub) vs multiplicity is shown in the 
right panel of Fig. [5] The measured V2 (~ 2.5%) in central U+U collisions is much lower than 
predicted (~ 4%) M35I1 . which requires further understanding and investigation of the production 
mechanisms for multiplicity and elliptic flow. 

The left panel of Fig. |6]shows 7 as a function of centrality for U+U collisions at 193 GeV, 



in comparison with Au+Au collisions at 200 GeV 112 lL 12411 . The opposite sign correlations in 
U+U are still higher than the same sign, with yss consistent with those in Au+Au and 70s 
slightly lower than those in Au+Au. To reduce the mutual background, we study 70s - 7ss, and 
multiply it by the number of participants, N pan , to compensate for the dilution effect [71]. The 
right panel of Fig.|6]shows the signal (70s - 7ss) '^Vpan vs v'2 for different centralities in 193 GeV 
U+U and 200 GeV Au+Au collisions. In both U+U and Au+Au, the signal roughly increases 
with V2, seemingly following the background trend described by Eq.(4). The central trigger in 
U+U is supposed to disentangle the background contribution from the signal, since the magnetic 
field will be greatly suppressed and the measured signal will be dominated by the V2-related 

6 



— 80 







-4- 



* 200 GeV Au+Au 
193 GeV U+U 

Opposite Sign:cc=-fS 

-a— 



5" 



60 



® Same Sign: a=p 

STAR preliminary 
. i .... i .... i .... i . 



40 20 

% Most Central 



o 

T— 

X 



^■200 

x 



o 



100 



200 GeV Au+Au 4^1 JF 

193 GeV U+U f% 



4 



★' 



STAR preliminary 



7 8 
V 2 (%) 



Figure 6: (Color online) Comparison between Au+Au collisions at 200 GeV (21, 24] and U+U collisions at 193 GeV. 
(Left) y as a function of centrality. (Right) (70s - Tss) X N pa n vs Vi. The error bars are statistical only. The open box 
represents the systematic uncertainty due to the tracking capability under the high luminosity in RHIC year 201 1 . 



background. As a result, in 0-1% most central U+U collisions the signal disappears as expected 
by the Chiral Magnetic Effect, while V2 is still ~ 2.5%. 



4. Summary 

The Chiral Magnetic Wave is a proposed signature of the Chiral Symmetry Restoration in the 
hot and dense nuclear matter created in heavy ion collisions at RHIC/LHC energies. From 200 
GeV to 19.6 GeV Au+Au collisions, the Vi(A^) difference between if~ and n + is consistentwith 
the qualitative expectations from the Chiral Magnetic Wave picture, and the slope parameter 
follows a centrality dependence qualitatively similar to the theoretical calculations of the CMW. 
On the other hand, UrQMD can not reproduce the data. Systematic investigations will be carried 
out to include different particle species and different V2 methods. 

In search of Local Parity Violation and Chiral Magnetic Effect, we studied the charge- 
dependent three-point correlator. The difference between the opposite sign and the same sign 
correlations is present in Au+Au, Cu+Cu, U+U and Pb+Pb collisions at RHIC and LHC, and 
remains almost unchanged up to 2.76 TeV and down to 11.5 GeV. The signal, yos - Tss, seems 
to disappear when the beam energy is lowered to 7.7 GeV, or when the magnetic field is greatly 
suppressed as in 0-1% most central U+U collisions, while V2 is still sizeable in both cases. The 
results seem to indicate that the observed signal is not dominated by the v'2-related background. 
A cross-check with 0-1% most central Au+Au collisions will be carried out in future. 
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